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Abstract 

Recently it was shown[l] that the imaginary part of the canonical parti- 
tion function of Schwarzschild black holes with an energy spectrum E n = 
a^/riEp ,n = 1,2, ... , has properties which - naively interpreted - leads to 
the expected unusual thermodynamical properties of such black holes (Hawk- 
ing temperature, Bekenstein-Hawking entropy etc.). 

The present paper interprets the same imaginary part in the framework of 
droplet nucleation theory in which the rate of transition from a metastable 
state to a stable one is proportional to the imaginary part of the canonical 
partition function. The conclusions concerning the emerging thermodynam- 
ics of black holes are essentially the same as before. 

The partition function for black holes with the above spectrum was calcu- 
lated exactly recently [1]. It is the same as that of the primitive Ising droplet 
model for nucleation in lst-order phase transitions in 2 dimensions. Thus 
one might learn about the quantum statistics of black holes by studying that 
Ising model, the exact complex free energy of which is presented here for 
negative magnetic fields, too. 
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1 Introduction 

In a recent paper I discussed properties of the canonical partition function 
for a quantum energy spectrum 



E n = a^/nE P ,n = 1,2,..., E P = c 2 y cK/G , a = O(l) , (1) 

where the n-th level has the degeneracy 

d n = <f , 5 > 1 . (2) 

The main interest in this spectrum comes from the many, differently justi- 
fied, proposals (see the corresponding quotations in ref. [1]) that a quantized 
Schwarzschild black hole might have such a spectrum! The associated canon- 
ical partition function 

oo , 

Z(t, x) = J2 e nt e~^ nx , t = \ng , x = (3aE P , (3) 

n=0 

converges only for \g\ < 1, whereas one is interested in its properties for 
g > 1. However, the function Z(g,x) can be continued analytically into the 
complex (yf-plane by means of the integral representation [1] 

x r°° e -z 2 /(4hiit) l 



Z(g = e\x) = —= / du — ^ (4) 

2^ Ji ln^w 



x r°° dr _ X 2/(4 T ) 



2v^7o r 3 / 2 i_ e (t-r) 



(5) 



which exhibits a branch cut of Z(<7, x) from g = 1 to g = +oo. Approaching 
the cut from above yields a complex Z, with an imaginary part 

and the principal value integral 

1 r 00 1 

Z r (t,x) =p.v.- drZi{T,x)- t — - (7) 

7T Jo 1 — e l ~ T 
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as its real part. 

Amazingly one gets essentially all the expected thermodynamical properties 
of a black hole, including the Hawking temperature, if one - naively - uses the 
imaginary part above for their derivation. This is strange and provoca- 
tive and needs further analysis. If such an approach can be substantiated, 
one will have an indication which quantum states of a black hole could be 
responsible for its unusual thermodynamical properties! 
In the following I shall try to put the properties mentioned into a certain per- 
spective and especially discuss their relationship to other approaches (which 
was not done in ref.[l]), namely the euclidean path integral (semiclassical) 
quantisation of black holes by Hawking and others and the theory of droplet 
nucleation associated with metastable states. 

The paper is organized as follows: In chapter 2 the properties of the model 
are briefly summerized and its theoretical and physical potentials and limita- 
tions sketched. Parts of it recall well-known material which serves to have the 
discussion more or less self-contained. Chapter 3 discusses the relationship 
of the model to Hawking's euclidean path integral approach to the thermo- 
dynamics of black holes, where he also gets a purely imaginary partition 
functional Chapter 4 deals with the nucleation of droplets occuring when 
a system changes from a metastable state to a stable one, e.g. in lst-order 
phase transitions or in the decay of "false" ground states. Langer|3[] was the 
first one to relate the imaginary part of the canonical partition function to 
the transition rate of such processes. Later Gross, Perry and Yaffe[Q in- 
terpreted the imaginary part of the semiclassical euclidean path integral - 
i.e. the canonical partition function - for the Schwarzschild black hole ("in- 
stanton") in this context, thus giving a different interpretation of the same 
quantity Hawking employed. 

In chapter 5 the relation of the partition function (3) to that of the primitive 
Ising droplet model for nucleation in 2 dimensions is discussed and it is very 
easy to see that both are exactly the same. The integral representation (4) 
therefore provides an exact free energy expression for this model, too, which 
as far as I know, is newQ. Chapter 6 presents some conclusions. 
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2 Scope of the model 



In the background of trying to push the quantum analysis - including its 
quantum statistics - of the simple isolated Schwarzschild gravitational system 
as far as possible are the following considerations: The quantized radiation 
emitted by a "classical" Schwarzschild black hole has a thermal (Planck) 
distribution governed by the (Hawking) temperature 

where M is the mass of the black hole. Assuming that this temperature 
of the radiation arises from its thermal contact with the black hole of the 
same temperature, the immediate question was - and still is -, how does the 
black hole gets this temperature? As the temperature (8) is proportional to 
Planck's constant the quantum theory of gravity has to play a decisive role in 
its understanding. Especially one would like to identify the microscopic quan- 
tum gravitational degrees of freedom which yield the corresponding quantum 
statistics and its macroscopic thermodynamics. 

Many attempts have been made to achieve this and I shall mention only a 
few of them in the following. The most active and enthusiastic reseach in 
this area presently takes place in connection with string theory. I shall say 
nothing about that in the following and refer to the very recent review by 
Horowitz 0. 

In view of the lack of a generally accepted quantum theory of Einstein's 
General Relativity it seems difficult to extract those microscopic gravita- 
tional degrees of freedom within this framework. However, the situation is 
not quite es hopeless as it might appear. The reason is that Einstein's the- 
ory for a rotationally symmetric isolated gravitational system, Schwarzschild 
gravity, can consistently be quantized, namely by first identifying its clas- 
sical observables in the sense of Dirac by solving the constraints associated 
with the gauge (diffeomorphims) degrees of freedom and then quantizing the 
remaining gauge-independent physical degrees of freedom afterwards. This 
has been done by Thiemann and myself|| in the framework of Ashtekar's 
formalism and briefly afterwards by Kuchaf|| in the geometrodynamical 
framework, the results being equivalent. They have been briefly summerized 



in ref.llO 



The essential point is that this system classically has just one canonical pair 
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of (Dirac) observables, namely its mass M and as the canonically conjugate 
quantity of M a time functional T[g; E] (of the metric g) which describes the 
difference r + — r_ in proper time r of two observers at the two (asymptotic) 
ends of the 1-dimensional spacelike hypersurface S which, together with the 
2-dimensional spheres S 2 , provides a time slicing of the 4-dimensional mani- 
fold. Formally the quantity 

* = Tfo; E] - (r+ - r_) (9) 
is the constant Dirac observable. 

Thus, for an observer at the asymptotic end where r —>■ +00 on the Schwarzschild 
manifold with the time r = r + there are - in the framework of this model - 
only two quantities available in order to describe the system: its Mass M and 
his proper time r! Everything else has to be expressed by them (and some 
fundamental constants like G, c, h, fee). Examples are: the Schwarzschild ra- 
dius Rs = 2MG/c 2 or any multiple thereof, the area A = AirR^ of the hori- 
zon, but also any time interval A associated with the system, i.e. A = 
where 7 is some number, e.g. 7 = Air for the time period A# of the euclidean 
section of the complex Schwarzschild manifold. 

Having eliminated the (infinite) gauge degrees of freedom classically one can 
now quantise the remaining physical degrees of freedom[8-10]. The extremely 
simple Schrodinger equation for this system, 

ihd T (j)(T) = Mc 2 (j)(T) , (10) 

has the plane wave solutions 

--Mc 2 r 

0(M,r)= X (M)e % , (11) 

where M is continuous and > 0. So at first sight there is no such spectrum 
like (1). However, one can get it by changing the boundary conditions. 
My doing so in ref. [10] was stimulated by the paper of Bekenstein and 
Mukhanov||ll||, in which they discussed properties of the spectrum (1) with 
a 2 = In (7/(471"), g = 2. Already in 1974 Bekenstein JTjJ] discussed this spec- 
trum as arising from a Bohr-Sommerfeld type quantisation of the area of the 
horizon (A oc n%) . Since then many authors have proposed such a spectrum 
(see refs. [3-23] in ref. [1]). 
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One knows from elementary quantum mechanics how to make the continuous 
momentum in plane waves discontinuous by periodic boundary conditions in 
space. Applying corresponding boundary conditions in time (with period A) 
to the plane waves (11) yields 

c 2 MA = 2nhn , n = 1, 2, . . . . (12) 

Introducing such boundary conditions is easier than justifying them by phys- 
ical arguments. In any case it means that the system is represented by the 
plane waves only for the time A after which it is abruptly terminated, e.g., by 
forming of the horizon. This may happen, e.g., because the state described 
by the plane wave is a metastable one - see below. There are several physcal 
time scales A = 'yRs/c associated with a black hole all of which[10] have 
7 ~ 0(1), so that 

E n = J-y/nE P = <Ty/nEp . (13) 



3 Euclidean path integral as canonical parti- 
tion function of a Schwarzschild black hole 

First I want to recall how Hawking's path integral approach to the canonical 
partition function in a semiclassical approximation leads to a purely imagi- 
nary partition function, too: 

Shortly after Hawking's discovery that the quanta emitted by a collapsing 
black hole have a thermal distribution governed by the temperature (8) Gib- 
bons and Hawking[13] employed the relationship between the canonical par- 
tition function of a system and the corresponding euclidean path integral in 
order to tackle the canonical quantum statistics of a black hole itself. For 
pure gravity that partition function Z is formally given by 

Z = f D\g]e- Ijs teV n , I E = —— I d 4 xR{g){g) 1/2 + surface terms, (14) 
J 16nG J 

where Ie is the euclidean action of pure gravity (for more details of the path 
integral approach to quantum gravity see the reprint collection in ref . [2] ) . In 
principle one has to take into account all field configurations g^ v which are 
periodic in euclidean time with period h/3 and which obey appropriate spatial 
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boundary conditions. In practice the path integral has been evaluated only 
semiclassically: One puts 

o 

sv =V +/ v > ( 15 ) 

o 

where 9^ is a classical solution of the euclidean field equations compatible 
with the boundary conditions just mentioned and the denote (small) 
quantum fluctuations around the classical solution. Inserting (15) into the 



action (14) and keeping only terms up to order 0{h 2 ^ u ) gives the semiclassical 
approximation 



if = + if , if = I d±x{g) 1/2 {hA{g)h) , (16) 



o 



where A{g) is a second order differential operator depending on the classical 
o 

metric 9^. 

Gibbons and Hawking evaluated Jg for the euclidean Schwarzschild metric 
and obtained 

Jf [Schwarzschild] = -(3 H Mc 2 , (17) 

where the period H/3h follows from the requirement that the euclidean Schwarzschild 
section has no conical singularities. The associated partition function 

Z ci = e -p H Mc 2 /2 = e -E 2 P p 2 H /(16n) (lg) 

gives the desired thermodynamics of black holes, especially the Bekenstein- 
Hawking entropy 

S B H/k B = A , £ = Gh/c 3 . (19) 

The approach appears to become "unconventional", if one includes the term 
if for which the path integral (14) is a Gaussian one which - formally - is 
proportional to the inverse square root of the product of the eigenvalues of 
the operator A (=(det(A)) -1 / 2 ). For the Schwarzschild case the term if 
has been evaluated by Gibbons and Perry[[14j]. A crucial point is that the 
operator A has just one negative eigenvalue[4] which contributes a factor i to 
the partition function. The combined contribution to the partition function 
Zg d of the classical euclidean Schwarzschild solution and the fluctuations 
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around it is given by M, WA 



Zg cl — iZ g Z G HP , 
ZcHP = l_^.( 0& )212/45 e -£M/(16,) iZa = e ,V/( C '^ 2p) 

where V is the spatial volume of the system and a is a "cutoff" associated 
with the ^-function regularization of det(v4). (The index "GHP" stands for 
"Gibbons, Hawking, Perry".) The factor Z g is the same as one would obtain 

for flat space (9/_ w = 5^). It represents the gas of free thermal gravitons 
surrounding the black hole. 

Except for the spatial volume factor the structure of Zghp in eq. (20) is very 
similar to that of Z{ in eq. (6) above: If a 2 = t/(4ir) - see ref. [1] - then the 
exponentials in both cases are equal for ft = ftu- The factors with the powers 
of ft - which in eq. (20) comes from short-distance fluctuations around the 
classical solution - are different in both cases. However this is not surprising 
because both approaches are so different. These factors may give rise to 
logarithmic corrections [1] to the Bekenstein-Hawking entropy (19) which is 
dominated by the exponential. I shall come back to this in the next chapter. 
Hawking has argued[2, ch. 15.8] that the partition function (20) has to be 
imaginary in order for the density of states N(E) in the Laplace transforms 

roo n T-i 1 /"-Moo 

Z(ft) = / dEN{E)e-f jE , N(E) = — / dftZ(ft)e? E (21) 
Jo 2ni J-ioo 

to be real if Z cl from eq. (18) is inserted into the last integral and which 
exists only if the contour of integration is rotated by 7r/2! Because of this 
Hawking has concluded || that the canonical ensemble is inappropriate for 
the description of black holes. 



Recently Hawking and Horowitz have discussed [15] a more refined version 
of the relationship between the classical euclidean gravitational action (14) 
- especially its boundary terms- and the entropy of the system. This work 



is related to that of York and collaborators |TB| who suggested to improve 
the approach of Gibbons and Hawking[13] and Gross, Perry and Yaffe[4] by 
employing a more microcanonical framework. Related is also the work of 
Iyer and Wald[|Hj . 



In any case, the imaginary parts Z GHP of eq. (20) or Zi of eq. (6), respec- 
tively, do seem to yield interesting though unconventional thermodynamical 
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properties, if one treates them as if they were "normal" partition functions. 
As they imply negative heat capacities they signal instabilities of the system. 
This leads us to a closely related interpretation of the imaginary part of a 
partition function, namely of being associated with the decay of a metastable 
state into stable one in the form of nucleation: 



4 Imaginary part of partition functions and 
droplet nucleation 

In a series of very influential papers Langer discussed the imaginary part 
of partition functions for systems with a lst-order phase transition where the 
transition from a metastable state to a stable one is initiated by the "nucle- 
ation" of droplets. Langer suggested that the nucleation transition rate is 
proportional to the imaginary part of the partition function. Coleman and 
Callen fil| discussed this approach in the context of instanton solutions in 



the euclidean version of otherwise lorentzean field theories and their role for 
the decay of "false" (metastable) vacua. The methods involved have found 
wide applications in different fields of physics |1| . 



The main procedure is very similar to the one described above for the eu- 
clidean section of the Schwarzschild solution: The instantons (or "bounces") 
are classical solutions of the " euclideanized" field equations satisfying ap- 
propriate boundary conditions of a certain path integral. The quantum (or 
thermal) fluctuations around these classical solutions leave the system stable 
if all the eigenvalues of the second- variation operator A (see above) are posi- 
tive. However, if one of the eigenvalues is negative, then the classical solution 
is not related to a minimum of the euclidean action integral, but only to a 
saddle point and the fluctuations can drive the system from its metastable 
state into a more stable lower one. 

The transition rate is essentially determined by ^s(Z scl ) /Zq, where Zq is the 
real part of the canonical partition function (which represents the metastable 
state) and by the absolute value of the negative eigenvalue. 
It is this framework in which Gross, Perry and Yaffe[f§] interpreted the contri- 
bution (20) of the euclidean Schwarzschild solution ("Schwarzschild instan- 
ton"): 

Eq. (20) represents the " 1-instanton" contribution to the partition function. 



9 



The next step is a dilute-gas approximation: Neglecting the "interactions" 
between the instantons, the contribution of N of them to the grand partition 
function Z G is 

Z g ±(iZ OHP F . (22) 

Summing over all N gives 

Z G = Z g e iZ GHP (23) 

and the grand canonical potential 

#09, V,a) = \nZ G = (3pV , = -Ud(3 + p(3dV - Nda , (24) 

(a = —fi(3, fi '■ chemical potential, p: pressure) here has the form 

* = \nZ g + iZ GHP , (25) 

whereas the general structure is [3, 19] 

* = lnZ + i%(Z sd )/Z . (26) 

Because of the special form of Z scl in (20) the factor Z g , representing the 
metastable graviton gas (i.e. Z g = Z ) drops out in eq. (26)! 
According to Langer|| (and others [18|, |T^]) the rate T of transition per unit 
volume from the metastable state (here the graviton gas) to the stable state 
(here the black hole) is given by 

T = = Mq.^/Zo) , V = */V , (27) 

7T 7T 

were k is the single negative eigenvalue mentioned above (that a large class of 
systems with possible negative eigenvalues has just one of them was proven 
by Coleman p0|1). 

The factor |«| in eq. (27) is a dynamical one, depending on nonequilibrium 
properties of the system[3,19]. As to more recent evaluations of \k\ for rela- 
tivists systems see refs. pTl . 

The relation between the entropy (19) and the exponential factor exp(— (3 2 / (lfwr)) 
in the transition rate (27) associated with eq. (25) is perhaps not so obvious 
anymore. However, in the case of nucleation one writes [3, 19] 

%{Z sd )/Z = e~^ pSd ~ F °) , F scl = -- \nQ{Z scl ) , F = -- lnZ (28) 

jj jj 
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and interprets F = F scl — Fq as the excess free energy of the critically large 
droplet. In our case we have Z = Z g and therefore Fs = —{1/(3) In Zqhp- In 
this sense Zqhp represents the partition function of the bare black hole, i.e. 
it describes the thermodynamics of the black hole without the surrounding 
graviton gas. 

Thus, the more sophisticated nucleation picture essentially leads to the same 
thermodynamical properties of Schwarzs child black holes as the "naive" use[l] 
of the imaginary part of the partition function! 

Still, it seems desirable to have a more systematic analysis of the relationship 
between the two approaches. 

Pagep2| has given a detailed analysis of the rate (27) for the transition of a 
gas into a black hole. The picture that black holes are formed by nucleation 
from a gas was already discussed by Gibbons and Perry [f23|j. 



5 Schwarzschild black hole quantum statis- 
tics and Ising droplet nucleation in 2 di- 
mensions 

One of the most amazing parts of the whole analysis presented here is that 
the quantum partition function (3) associated with a Schwarzschild black 
hole is the same as that of the (primitive) Ising droplet model for nucleation 
(in lst-order phase transitions) in 2 dimensions ^ ]: 

Suppose a d-dimensional (d > 2) lattice of N Ising spins below the critical 
temperature T c to be in a positive magnetic field H so that almost all the 
spins are up (+1). If one then slowly turns the magnetic field negative, the 
system comes into a metastable state in which the total magnetization is still 
positive, before finally becoming negative, too. A simple model to describe 
the transition into that stable state is the following: Assume that (small) 
droplets containing / down-spins form, in a background of up-spins. The 
droplets are supposed to be noninteracting (dilute-gas approximation) and 
the number n\ of droplets of size I to be given by 

m = Ne~P €l ,1 = 1,... . (29) 
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The droplet formation energy e; is assumed to consist of two (competing) 
terms, the bulk energy 2HI and a surface energy term (f)l^ d ~ 1 >' d J so that 

e t = 2HI + <j)l( d ~ l ">/ d . (30) 
The interesting part of the canonical partition function is the finite sum 

Z N /N = J2e-^ Hl -^ l{d ~ 1)/d . (31) 
i 

Letting N go to oo and employing the grand canonical ensemble in exactly the 
same way as for the instanton gas above yields the following grand canonical 
potential if) per spin 



mH)=J2e- 2lH Pe- lid - 1)/d ^ . (32) 

1=0 

The theoretical investigations of the function if) have focussed mainly on 
its behaviour for negative H (where the sum no longer converges!) in the 
neighbourhood of H = 0. It is clear that if) must have a singularity at 
H = 0, the so-called "condensation point". Many approximate analytical 
calculations and rigorous estimates ||25|| suggested that if) has an essential 
singularity there, a supposition supported by numerical calculations ||26|| . 
Now for d = 2 the series (32) is just the same as (3) which has been summed 
exactly in ref. [1] by using Lerch's observation^] that the relation 



\x\ 


POO 




1 dve~ 


V^J 





\x\ 


r°° dr 


2^ 


Jo r 3 / 2 



_V^2 = \x\_ [°° dvft -a? v 2 /4:-n/v 2 (33) 

-x 2 / (At) — tit 

turns the series into a geometrical one which can be summed under the 
integral sign and then continued analytically. We only have to put t = 
—2/3H, x = (f>(3 and then get from eq. (6) the exact result 

S(^)(/3,iY< 0) = ^^^/5- 1 /2| iT |-3/2 e -(^ 2 /?)/(8| J f/|) (34) 
4v 2 

Thus, there is indeed an essential singularity at H = 0. Instead of |i/|~ 3 / 2 
in front of the exponential factor the approximation methods of Langer[3] 
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and Giinther et al.[25] yield the factor \H\. The modified method used by 
Harris [26] gives the correct \H\ -dependence! 

The real part of ip is to be calculated by means of the principal value (dis- 
persion) integral (7) (if the sum (32) starts with I = 1 etc. corresponding 
changes have to be madefl]). 

In ref.[l] I pointed out that the real part (7) can become negative for small 
x. The droplet nucleation picture might help to understand the (physical?) 
background of this. 

For a recent discussion of metastability in the 2-dimensional Ising model it- 
self see ref . . 

It was already stressed in ref.[l] that the series (3) obeys the heat equation 
dtZ = d 2 Z. The corresponding partial differential equation for if> in d=3 
dimensions is d 2 ip = —d^pp. This ought to help analysing the behaviour 
of ip in the neighbourhood of t = for d=3, too! I shall only indicate 
what to expect: As that partial differential equation is invariant under the 
scale transformation t — > At, x — > \ 2 ^ 3 x, there should exist solutions of the 
type f{y),y = x 3 /t 2 . Inserting this gives for f(y) the ordinary diff. eq. 
27 f" + (4 + 5A/y)f" + 6(1/?/ + l/y 2 )f = 0. For t -> 0, y -> oo, there is 
the approximate solution f" ~ exp(— (4/27)?/) which agrees with the form of 
the expected essential singularity as to small negative if in 3 dimensions [25]. 
However, more analysis is certainly necessary here. 



6 Conclusions 

It appears that there can be hardly any doubt that the square root spec- 
trum (1) is closely related to the expected thermodynamics of black holes. 
It is therefore important to understand its physical meaning and the theo- 
retical background better. Because it implies that the area A of the horizon 
is proportional to h one should probably look here for an interpretation [11]. 
This corresponds to other recent investigations which stress the importance 
of bifurcate horizons and the associated surfaces for understanding the ther- 
modynamics of black holes [[29|j . 

Furthermore, the discussion above shows that a "naive" and a droplet nucle- 
ation interpretation, respectively, of the imaginary part of a partition func- 
tion provide just two different aspects of the thermodynamical properties of 
metastable systems. 
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Similarly important is the physical understanding of the degeneracy (2) which 
is equally essential for the derivation of the thermodynamics. In the Ising 
droplet model above g is replaced by the factor exp —f3H which represents 
the influence of the driving external (negative) magnetic field. The same role 
can have a (positive) chemical potential [3, 19]. Bekenstein and Mukhanov[ll] 
used information theoretical arguments to put g = 2. However, the argu- 
ments above go through with any g > 1. The physical background of the 
degeneracies (2) probably needs further understanding. (The "driving field" 
behind g > 1 is most likely the gravitational attraction which leads to the 
black hole.) 

In addition it will be interesting to see whether and how the above consid- 
erations can be extended to the thermodynamics of the Reissner-Nordstr0m 



model |30] where investigations similar to those of refs. [8,9] for the Schwarzschild 



case exist [RT 



Very exciting is the possibility to map the quantum statistics of a Schwarzschild 
black hole onto the statistics of the 2-dimenionsal Ising model for droplet 
nucleation. This might allow to learn from a known field of physics for an 
unknown one. 

Finally it is surprising that the pursuit of black hole physics leads to an exact 
mathematical solution for an old condensed-matter problem. 

I am very much indebted to Malcolm Perry for drawing my attention to 
refs. [3] and [4]. 

Note added in response to a question of the referee: 



If one considers a Schwarzschild black hole in d+1 space-time dimensions 132 
d > 3, then its Schwarzschild radius Rs is proportional to M l / Iyd ~ 2 \ where 
M is the mass of the system. Thus, if we again assume the time interval A 
in eq. (12) above to be proportional to Rs, then ^/n in eq. (13) is replaced 
by n ( d_2 )/( d_1 ) and the resulting partition function is the same as that of the 
Ising droplet model for nucleation in d-1 dimenions (see eq. (32))! 
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